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The major goal of the various heavy ion programs is the search for a transient state 
of deconfmed matter, dubbed the quark-gluon-plasma (QGP): A phase transition to this 
new state of matter is predicted by lattice QCD when a sufficiently high energy density 
£^ ■ (e ~ 1 GeV/fm 3 ) is reached. 

Strange particle yields and spectra are key probes to study excited nuclear matter and 
to detect the transition of (confined) hadronic matter to quark-gluon-matter (QGP)[[l|. [|. 
The relative enhancement of strange and multi-strange hadrons, as well as hadron ratios 
in central heavy ion collisions with respect to peripheral or proton induced interactions 
have been suggested as a signature for the transient existence of a QGP-phase [ 

A wealth of systematic information has been gathered to study the energy dependence 
of observables from yfs ~ 2 AGeV to ^/s = 200 AGeV. For the first time also information 
on unstable particle emission in AA reactions is available: the $ and A(1520) have been 
^ ■ observed in heavy ion reactions at SPS energies [ [3J, SPS and RHIC experiments [ [J[ || 
report measurements of the if (892) signal, and RHIC is already attacking the fo and p 
mesons and the £(1385). 

A key problem is the identification of unambiguous signatures of possible QGP cre- 
ation. Under the assumption of thermal and chemical equilibrium, fits with a statistical 
(thermal) model have been used to extract bulk properties of hot and dense matter, e.g. 
the temperature and chemical potential at which chemical freeze-out occurs [ [5], || 0. In 
addition also dynamical evaporation models and non-equilibrium transport calculations 
have been employed to study the energy and centrality dependence of particle and espe- 
cially strangeness productionf ||, H, [K| [11]]. Overall, the conclusions are ambigous and 



range from evidence for QGP creation to canonical enhancement and from string fusion 
to baryon junctions. 

Numerous ideas have been brought forward in the past, ranging from the search for 
a softest point in the equation of state, photon and lepton radiation, J/\l/ suppression 
to event-by-event fluctuations. However, the main difficulty in the interpretation of the 
available data is that the observed final state hadrons carry relatively little information 
about their primordial sources. Most of the hadrons had been subject to many secondary 
interactions and are strongly influenced by the decays of high mass resonances. To shed 
some light on the strangeness production mechanism and the question of chemical and 
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kinetic equilibration and decoupling, we study 

• the production of multiple strange baryons in pp interactions. Here on can directly 
probe the microscopic decay of color flux tubes, allowing to differentiate between 
different string models and a statistical description of the hadronization. 

• The energy and centrality dependence of (strange) hadron resonances in AA which 
carry unique information about the stage between chemical and kinetic freeze-out. 

Hadronization in proton-proton interactions 

In the string picture, high energy proton-proton collisions create "excitations" in form of 
strings, being one dimensional objects which decay into hadrons according to longitudinal 
phase space. This framework is well confirmed in low energy electron-positron annihilation 
P~2|] where the virtual photon decays into a quark-anti-quark string which breaks up into 



various kinds of hadrons. However, specific string models differ in their philosophy and 
the types of strings that are created. In general two classes of models - based on color 
or momentum exchange - can be distinguished, the resulting objects two quark-diquark 
strings with valence quarks being their ends, however, is quite similar. Here we will 
contrast the present UrQMD prescription for string formation (similar to the PYTHIA 
model [ |H| with the approach used in NeXuS 3.0: In UrQMD [ [L4| the projectile and 
target protons become excited objects due to the momentum transfer in the interaction. 
The resulting strings, with at most two strings being formed, are of the di- quark- quark 
type. 
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Figure 1. a) The simpliest collision configuration has two remnants and one cut Pomeron 
represented by two q — q strings, b) One of the q string-ends can be replaced by a qq string-end. 
c) With the same probability, one of the q string-ends can be replaced by a qq string-end (taken 
from [|§). 



In NeXuS 3.0[ [T5|] , the pp interaction is described in terms of Pomeron exchanges or 
ladder diagrams. Both hard and soft interactions take place in parallel. Energy is shared 
equally between all cut Pomerons and the remnants. Here, for the string ends of soft and 
semi-hard Pomerons quarks and anti-quarks from the sea are taken in a flavour-symmetric 
way. Thus, the valence quarks stay in the remnants, yielding excited quark bags [ [18| (see 
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Fig. [I]). In contrast to these string models, the predictions of two statistical models (model 



I, being fully canonical [ [16] and model II, being canonical with respect to strangeness [ 
|17|1 ) are also shown. 

Fig. |2| (left) depicts the anti-baryon to baryon ratio at midrapidity in proton proton 
interactions at 160 GeV. The results of the calculations by the new NeXuS 3.0 (ratios cal- 
culated from [ and UrQMD/PYTHIA, which are well established string-fragmentation 
models for elementary hadron hadron interactions, are included in this figure. In both 
di-quark models, the B/B ratio increases strongly with the strangeness content of the 
baryon. For strangeness \s\ =3 the ratio significantly exceeds unity. In UrQMD and 
PYTHIA the hadronization of the di-quark- quark strings leads directly to the overpopu- 
lation of Q as discussed in detail in [ |l9fl . The new string formation scheme employed in 
NeXuS 3.0, however, allows to get a reasonable agreement with experimental data. The 
basic features of the production of multi-strange baryons as well as of As and protons 
can be understood within the model picture of proton-proton scattering: the created fi- 
nal particles emerge from a non-trivial system of projectile/target remnant states and a 
number of cut Pomerons each represented by a pair of strings. 
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Figure 2. Left: anti-baryon to baryon ratio at \y — y cm \ < 1 in pp interactions at 160 GeV 
as given by PYTHIA, NeXuS3.0 and UrQMD. Right: anti baryon to baryon ratio for the 
same reaction as given by statistical models. Stars depict preliminary NA49 data for the 
B/B ratio at midrapidity. 



The predictions of the statistical models (in full phase space) are shown in Fig. ^| (right). 
In these approaches the B/B ratio is seen to exhibit a significantly weaker increase with 
the strangeness content of the baryon than expected in the di-quark string fragmentation 
models. In the grand canonical picture, where the B/B ratio is very sensitive to the 
baryon chemical potential fiB, it is easily understood that statistical models, can not 
yield a ratio of Q/Q > 1. For finite baryon densities as in the pp system, the B/B has 
to stay below one and only in the limit of \is = may Q/Q = 1 be approached. This 
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Figure 3. Top: Inelastic and (pseudo-)elastic collision rates in Pb+Pb at 160AGeV. r c h 
and Tth denote the chemical and thermal/kinetic freeze-out as given by the microscopic 
reaction dynamics of UrQMD. Bottom: Average energy per particle at midrapidity (\y — 
Vcm\ < 0.1) as a function of time (taken from [ |23|). 



feature is not blurred in the canonical approach. For comparison, both figures include 
preliminary data on the B/B ratios obtained at midrapidity by the NA49 Collaboration 
[ (Very preliminary NA49 data seems to support a Q/Q ratio below one [ |21J). 

It is important to note that the result Q/Q > 1 in pp collisions from the di-quark string 
models solely depends on the geometry of the decaying objects. The ratio can not be 
altered (or even inverted) by modifications of the strangeness and di-quark suppression 
factors. These factors can only modify the absolute yields of strange hadrons, but do not 
influence the ratio discussed here. 

If the new NA49 data can be confirmed, one is forced to conclude that the di-quark 
string models fail to describe the ratios of multiple strange baryons in pp interactions at 
the SPS. A solution might be to replace the conventional approach by a system of sea 
quark strings and remnants or to abandon the flux tube picture, which has explained 
many dynamical features of hh collisions, for a statistical hadronization model. 

Chemical and kinetic freeze-out in AA collisions 

In nucleus-nucleus collisions the situation is by far more involved than in the pp case 
discussed above. Here, the major fraction of finally observed particles in the experimental 
setup stem from decays of resonances (mesonic or baryonic) which have undergone many 
scatterings from their point of production to observation. The final hadron yields seem 
to be compatible with a hadronic gas described by the baryo-chemical potential \ib and 
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a temperature parameter T in a statistical model (see e.g. [ H El El)- ^ nas been 
suggested that both parameters are coupled by a universal freeze-out criterion assuming 
a mean energy per hadron (E)/(N) = 1 GeV/hadron [ 2l| at chemical freeze-out. In 
this framework, the formed hadrons do only undergo elastic collision from this chemical 
freeze-out time (at SPS energies at a temperature of ~ 160 — 170 MeV to the final kinetic 
freeze-out at temperatures of the order of 120 MeV. 

To explore whether this sequential freeze-out is realized in heavy ion reactions at highest 
energies and how it can be probed, the Ultra-relativistic Quantum Molecular Dynamics 
model (UrQMD 1.2) is used [ [14]]. This microscopic transport approach is based on the 
covariant propagation of constituent quarks and di-quarks accompanied by mesonic and 
baryonic degrees of freedom. As discussed above, the leading hadrons of the fragmenting 
strings contain the valence-quarks of the original excited hadron and represent a simplified 
picture of the leading (di) quarks of the fragmenting string. The elementary hadronic 
interactions are modelled according to measured cross sections and angular distributions. 
If the cross sections are not experimentally known, detailed balance is employed in the 
energy range of resonances. The partial and total decay widths are taken from the Particle 
Data Group. 

To analyse the different stages of a heavy ion collision, Fig. |3| (top) depicts the time 
evolution of the elastic and inelastic collision rates in central Pb+Pb interactions at 
160AGeV. The inelastic collision rate (full line) is defined as the number of collisions with 
flavour changing processes (e.g. nn — > KK). The elastic collision rate (dashed dotted 
line) consists of two components, true elastic processes (e.g. Kn — > Kn) and pseudo-elastic 
processes (e.g. Kir — > K* — > Kit). While elastic collision do not change flavour, pseudo- 
elastic collisions are different. Here, the ingoing hadrons are destroyed and a resonance 
is formed. If this resonance decays later into the same flavours as its parent hadrons, 
this scattering is pseudo-elastic. Fig. |3| (bottom) depicts the average thermal energy - 
calculated from interacted hadrons with pi = [p 2 x +Py) /2 - per particle at midrapidity 

(\y-y cia \ <0.l). 

The present microscopic study of the collision dynamics, indeed supports the idea of 
separated phases in the evolution of the system: 



I t < 2 fm/c: In the initial stage of the nucleus-nucleus reaction non-equilibrium 
dynamics leads to strong baryon stopping in multiple inelastic interactions, shown 
by huge and strongly time dependent collision rates. This stage deposits a large 
amount of (non-thermalized) energy and creates the first generation particles. 

II 2 fm/c< t < 6 fm/c: Due to the high particle densities and energies inelastic 
scatterings dominate this stage of the reaction. Chemical equilibrium might be 
achieved due to a large number of flavour and hadro-chemistry changing processes 
until chemical freeze-out. 



Ill 6 fm/c< t < 11 fm/c: After the system has expanded and cooled down, elastic and 
pseudo-elastic collisions take over. Here, only the momenta of the hadrons change, 
but the chemistry of the system is mainly unaltered, leading to the kinetic freeze-out 
of the system. 
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IV t > 11 fm/c: Finally, the reactions cease and the scattering rates drop drastically. 
The systems breaks up. 

Even the hypothesis of chemical decoupling at an energy per hadron of 1 GeV, is in 
line with our analysis as demonstrated in Fig. |3| (bottom). In the non-equilibrium stage 
(E) I (N) decreases steadily. While the system is 'cooking' in the inelastic scattering stage, 
(E) I (N) stays constant around 1 GeV and drops suddenly as the hadronic system enters 
the kinetic stage. 

The spectra and abundances of A(1520), K°(892) and other resonances can be used to 
study the break-up dynamics of the source between chemical and thermal freeze-out. If 
chemical and thermal freeze-out are not separated - e.g. due to an explosive break-up of 
the source - all initially produced resonances are reconstructable by invariant mass analysis 
of the final state hadrons. However, if there is a separation between the different freeze- 
outs, a part of the resonance daughters rescatter, making this resonance unobservable in 
the final state. Thus, the relative suppression of resonances in the final state compared 
to those expected from thermal estimates provides a chronometer for the time period 
between the different reaction stages. Even the chemical composition of the system might 
be changed by up to 10% in the hyperon sector (after chemical 'freeze-out'), due to 
inelastic scatterings of the resonance daughters (e.g. Kp — ► A). 

The rapidity spectra for A(1232), A*(1520), K°(892) and $ in Pb(160 AGeV)Pb, b< 
3.4 fm collisions, are depicted in Fig. |j. The left part, shows the total amount of decaying 
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Figure 4. Rapidity densities for A(1232), A*(1520), K°(892) and $ in Pb(160 AGeV)Pb, 
b< 3.4 fm collisions. Left: All resonances as they decay. Right: Reconstructable reso- 
nances (taken from [ |23| ) . 

resonances. Here, subsequent collisions of the decay products have not been taken into 
account - i.e. whenever a resonance decays during the systems evolution it is counted. 
However, the additional interaction of the daughter hadrons disturbs the signal of the 




7 



resonance in the invariant mass spectra. This lowers the observable yield of resonances 
drastically as compared to the primordial yields at chemical freeze-out. Fig. ^ (right) 
addresses this in the rapidity distribution of those resonances whose decay products do not 
suffer subsequent collisions - these resonances are in principle reconstructable from their 
decay products. Note that reconstructable in this context still assumes reconstruction of 
all decay channels, including many body decays. 

By using the estimates done by [ Q in a statistical model, it seems possible to relate the 
result of the present microscopic transport calculation to thermal freeze-out parameters. 
The surprising result is that the microscopic source seems to have a lifetime shorter than 
1 fm/c and a freeze-out temperature lower than 100 MeV. Apparently, the values obtained 
from UrQMD seem to favour a scenario of a sudden break-up of the initial hadron source, 
in contrast to the time evolution of the chemical and thermal decoupling as shown in Fig. 
[3[ This misleading interpretation might be due to the re-creation of resonances in the 
elastic scattering stage, which was not taken into account in the statistical model analysis. 
The influence of these effects is currently under investigation. 

However, even the question of the existence of such resonance states in the hot and 
dense environment is still not unambiguously answered. Since hyperon resonances are 
expected to dissolve at high energy densities (see .e.g. [ |25|) it is of utmost importance 
to study the cross section of hyperon resonances as a thermometer of the collision. 

To set the stage, Figs. [| and |6] address the excitation function of observable resonance 
multiplicities. In addition, a comparison of rapidity integrated yields (4-7T values, circles) 
with the center-of-mass values (y cm ±0.5, squares) is given. The A includes decays from 
S°, but not from S. Protons do not include decays from A's. All strong decays are 
included. No cuts are applied except when mentioned. The anti-K* multiplicities are 
monotonously increasing with energy, while the hyperon resonance show a pronounced 
maximum in the excitation function at y/s =8 GeV (Ei a b ~ 30 AGeV). This maximum 
is also present - as earlier observed for hyperons by [ |26| - if scaled by the number of 
pions. This makes the newly planned GSI200 facility an ideal place to study in-medium 
modifications of resonances. 

More information can be obtained if the multiplicities are normalised to the groundstate 
hadrons, i.e. (anti-)Kaons and Lambdas. Figs. [7| and || show the energy dependence of 
the h*/h ratios. Here the baryo-chemical potentials cancel out and information on the 
freeze-out temperature can be obtained. 

To summarise, hadronization into multiple strange hadrons in elementary pp inter- 
actions has been studied within different string approaches and statistical models and 
confronted with data. Di-quark string models seem not to be in line with the experimen- 
tal observations of Q and Q production in pp at the SPS. The freeze-out dynamics in 
AA has been microscopically explored. Different decoupling stages have been identified 
in the non-equilibrium dynamics. Observables to clock the kinetic scattering stage and 
the freeze-out temperature with (strange-) resonances were discussed. 
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Figure 5. Multiplicity excitation function for central Pb+Pb (Au+Au) reactions, in 4n 
(circles) and at midrapidity (squares). Left: (if Q (892)). Right: (Z°(892)). 
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